PROPAGATION OF CORRELATIONS 



IN QUANTUM LATTICE SYSTEMS 

BRUNO NACHTERGAELE, YOSHIKO OGATA, AND ROBERT SIMS 

Abstract. We provide a simple proof of the Lieb- Robinson bound and use it to prove the existence 
of the dynamics for interactions with polynomial decay. We then use our results to demonstrate 
that there is an upper bound on the rate at which correlations between observables with separated 
support can accumulate as a consequence of the dynamics. 



1. Introduction 

Recently, there has been increasing interest in understanding correlations in quantum lattice 
systems prompted by applications in quantum information theory and computation ^3 [3 03 IH and 
the study of complex networks jK] • The questions that arise in the context of quantum information 
and computation are sufficiently close to typical problems in statistical mechanics that the methods 
developed in one framework are often relevant in the other. The bound on the group velocity in 
quantum spin dynamics generated by a short-range Hamiltonian, which was proved by Lieb and 
Robinson more than three decades ago [S], is a case in point. For example, as explained in j2], the 
Lieb-Robinson bound provides an upper bound on the speed of information transmission through 
channels modeled by a quantum lattice systems with short-range interactions. 

The Lieb-Robinson bound plays a crucial role in the derivation of several recent results. For some 
of these results it was useful, indeed necessary, to generalize and sharpen these bounds. Several 
such improvements have recently appeared [§1 El- in this paper we provide a new proof of the 
Lieb-Robinson bound (Theorem 12. 1|) and other estimates based on a norm-preserving property of 
the dynamics (see Lemma |A.1|) , We apply this result to give upper bounds on the rate at which 
correlations can be established between two separated regions in the lattice for a general class of 
models (Theorem K-i.ljl , Moreover, our bounds allow us to prove the existence of the dynamics 
(Theorem 12. 2|) , in the sense of a strongly continuous group of automorphisms on the algebra of 
quasi-local observables for a larger class of interactions than was previously known [T| 17]. 

1.1. The Set Up. We will be considering quantum spins systems defined over a set of vertices A 
equipped with a metric d. A finite dimensional Hilbert space 7i x is assigned to each vertex x £ A. 
In the most common cases A is a graph, and the metric is given by the graph distance, d(x,y), 
which may be the length of the shortest path of edges connecting x and y in the graph. 

For any finite subset X C A, the Hilbert space associated with X is the tensor product Tlx = 
®xex T^x, and the set of corresponding observables supported in X is denoted by Ax = B(Ttx), the 
bounded linear operators over Tlx- These local observables form an algebra, and with the natural 
embedding of Ax x in Ax 2 for any X\ C X2, one can define the C*-algebra of all observables, A, as 
the norm completion of the union of all local observable algebras Ax for finite X C A. 

An interaction is a map <£ from the set of subsets of A to A with the property that &(X) € Ax and 
&(X) = &(X)* for all finite X C A. A quantum spin model is then defined to be the Hamiltonian, 
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expressed in terms of its interaction, given by 

(1.1) H* := *(*)• 

XcA 

For notational convenience, we will often drop the dependence of H$ on $. 

The dynamics, or time evolution, of a quantum spin model is the one-parameter group of auto- 
morphisms, {Tt}i g iR, defined by 

(1.2) n{A) = e ltH Ae- itH , A £ A, 

which is always well defined for finite sets A. In the context of infinite systems, a boundedness 
condition on the interaction is required in order for the finite-volume dynamics to converge to a 
strongly continuous one-parameter group of automorphisms on A. 

To describe the interactions we wish to consider in this article, we first put a condition on the set 
A; which is only relevant in the event that A is infinite. We assume that there exists a non-increasing 
function F : [0, oo) — > (0, oo) for which: 
i) F is uniformly integrable over A, i.e., 

(1.3) ||F|| := sup^F(d(x,y)) < oo, 



and 

ii) F satisfies 



xeA yeA 



^F(d(x,z)) F(d(z,y)) 



(1.4) C := sup > v v ' /;. \; < oo. 

Given a set A equipped with a metric d, it is easy to see that if F satisfies i) and ii) above, then 
for any a > the function 

(1.5) F a (x):=e- ax F(x), 

also satisfies i) and ii) with ||-F a |i < H-^ll an d C a < C. 

As a concrete example, take A = Z rf and d(x, y) = \x — y\. In this case, one may take the function 
F(x) = (1 + x)~ d ~ e for any e > 0. Clearly, (|1.3|) is satisfied, and a short calculation demonstrates 
that (jOjl holds with 

(!-6) C < 2 d+£+1 Y . 

We also observe that, although the purely exponential function G{x) = e~ ax , is integrable for a > 0, 
i.e., it satisfies i), it does not satisfy ii). This is evident from the fact that the cardinality of the 
set {z € 1j d : \x — z\ + \z — y\ — \x — y\ = 0} is proportional to \x — y\, and therefore, there exists 
no constant C uniform in \x — y\. 

To any set A for which there exists a function F satisfying i) and ii) above, we define the set 
B a (A) to be those interactions $ on A which satisfy 

(i-7) \m\ a ■■= sup JJ}?\ < °°- 

2. Lieb- Robinson Estimates and Existence the Dynamics 

2.1. Lieb-Robinson Bounds. We first present a variant of the Lieb-Robinson result which was 
first proven in [<3 El • 
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Theorem 2.1 (Lieb-Robinson Bound). Let a > and take A\ C A a finite subset. Denote by T t 1 
the time evolution corresponding to a Hamiltonian 



(2.1) 



XcAi 



defined in terms of an interaction $ £ B a (A). There exists a function g : R — > [0, oo) wii/i t/ie 
property that, given any pair of local observable A 6 .Ax and 1? 6 Ay with X,Y C Ai, one may 
estimate 

2\\A\\ \\B\ 



(2.2) [r^{A),B] 

for any t € R. .Here t/ie function 
(2-3) s a (t) = 



< 



( e 2||*||„C |t|_ 1 ) ifd(X,Y)>0, 
e 2 H*H a Ca I*' otherwise. 



Proof. Consider the function / : R — > A defined by 

(2.4) /(*):= [7^(i4),J3]. 
Clearly, / satisfies the following differential equation 

(2.5) f'(t) = i [f(t),r t A i (H X )] +i [t^(A), [r t M (H x ),B 

where we have used the notation 

(2.6) Hy= 

ZCA 1 : 

for any subset Y C A\. The first term in ()2.5|) above is norm-preserving, and therefore the inequality 



(2.7) 



[rt(A),B}\\ < \\[A,B]\\+2\\A\\ I *' \\[r^(H x ), B]\\ ds 

Jo 



follows immediately from Lemma I A. II and the automorphism property of r f Al . If we further define 
the quantity 

(2.8) C bP M) := sup i^'^ 

A£Ax 



\A\\ 



then ()2.7() implies that 
(2.9) C B (X,t) <C B (X,0)+2 V \MZ)\\ C B (Z,s)ds. 



ZCA X : 



Clearly, one has that 
(2.10) 



C B (Z,0) < 2\\B\\5 Y (Z), 



where 5y{Z) = if Z C\Y = and S Y (Z) = 1 otherwise. Using this fact, one may iterate Q2.9|) and 
find that 



(2.11) 



C B (X,t) < 2\\B\\ 



(2\t\ 



n=0 
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where 

n 

(2-12) a n = E E 

Z 1 CA 1 : Z 2 CA 1 : Z„CA 1: i=l 

Zinx^0z 2 nZi^0 z n nz n _i^0 
For an interaction * G £> a (A), one may estimate that 

(2.i3) °i < EE E 11^)11 ^ n*iu EE f »^»))" 

In addition, 

«2 < EE E E u*( z i)ii E n $ (^)H 

iSXl/gy ZgAl ZjCA^ Z 2 CA i: 
Zi3z,z Z2Bz,y 

< wnl EEE F «(^^))^(^y)) 

x£X yGY zGA 

(2-14) < C a ££F a (d(x,y)), 

xGX j/GY 

using (|1.4j) . With analogous arguments, one finds that 

(2.15) a n < ii*ii™ cr 1 EE^^-f))' 

x£X y£Y 

Inserting (|2.15j) into (|2,11|) we see that 

(2.16) C B (X,t) < ^M exp [2||cI>|| a C a !t|] ^2^F a (d(x,v)), 

a x£XyeY 

from which (|2.2|l immediately follows. 

In the event that d(X, Y) > 0, one has that Cb(X, 0) = 0. For this reason the term corresponding 
to ao = 0, and therefore, the bound derived in ()2.16|) above holds with e 2 "*"™^'*' replaced by 

e 2||*|| a C„|t| _ 1. □ 

We note that, for fixed local observables A and B, the bounds above are independent of the 
volume Ai C A. 

In the event that * € B a (A) for some a > 0, then the bound in Q2.2j) implies that 
(2-17) [r^(A),B] < ^^||F||min(|X|,|y| )e -»[^-^^, 

which corresponds to a velocity of propagation given by 

2||*LG I 



(2.18) V$ := inf 

a>0 a 

We further note that the bounds in (|2.2j) and (|2.17j) above only require that one of the observables 
have finite support; in particular, if \X\ < oo and d(X, Y) > 0, then the bounds are valid irrespective 
of the support of B. 

One can also view the Lieb-Robinson bound as a means of localizing the dynamics. Let A be 
finite and take X C A. Denote by X c = A \ X. For any observable A 6 Aa set 



(2.19) {A) x *:= U*AUfi{dU), 

JU(X C ) 
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where U(X C ) denotes the group of unitary operators over the Hilbert space TLx c and \i is the 
associated normalized Haar measure. It is easy to see that for any A E Aa, the quantity {A)x c E Ax 
and the difference 



(2.20) 



(A) 



X'~ 



A 



U{X C ) 



U* [A,U] nidU). 



We can now combine these observations with the Lieb-Robinson bounds we have proven. Let 
A E Ax be a local observable, and choose e > 0, a > 0, and an interaction <& E B a (A). We will 
denote by 

B t {e) = B(A,t,e) := \ x E A : d(x,X) < 2 ^ a ° a \ t \ + e 



(2.21) 

the ball centered at X with radius as specified above. For any U E U(Bf(e)), we clearly have that 
(2.22) 



d(X,supp(C/)) > 2 ^ aCa \ t \ + e , 



and therefore, using 1)2.20(1 above, we immediately conclude that 

\n(A) - {n(A)) B}{£) 

(2.23) < 
where for the final estimate we used (|2,17|) . 



< / \\[r t (A),U]\\ n{dU) 

U{B1{e)) 
2\\A\\ \X\ 



C a 



File 



2.2. Existence of the Dynamics. As is demonstrated in pQ, one can use a Lieb-Robinson bound 
to establish the existence of the dynamics for interactions <3? E B a (A). In the following we consider 
the thermodynamic limit over a increasing exhausting sequence of finite subsets A n C A. 

Theorem 2.2. Let a > 0, and <3? E B a (A). The dynamics {rt}t£R corresponding to $ exists as a 
strongly continuous, one-parameter group of automorphisms on A. In particular, 



(2.24) 



lim ||r t A "(A)-Tt(^)|| =0 



for all A E A. The convergence is uniform for t in compact sets and independent of the choice of 
exhausting sequence {A n }. 

Proof. Let n > m. Then, A m C A n . It is easy to verify that for any local observable A E Ay, 

rt d 



(2.25) 
and therefore 
(2.26) 



r^{A)-r^{A) = J ± (r A "r^(A)) ds, 

t*»(a) -r t A -(A)\\ < y, e r ii i*w> r * Am w ii ds - 

_ A \ A XT —«_-*' 



x£A n \A m Xlx ' 

Applying Theorem 12. H we see that the right hand side of (|2.26|) is bounded from above by 

-1*1 



(2.27) 



2|L4|| f g a (s)ds Yl EH $ ( X )HEE^(^y))- 



x£A n \A m XBx 
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Rewriting the sum on X 3 x and y & X as the sum on y G A and X B x,y, one finds that 

r\t\ 

t^(A)-t^(A) < 2|L4||||$|| a C a / g a (s)ds £ £ F a (d(x, z)) 

xeA„\A m zeY 

1*1 



(2.28) < 2\\A\\\\<S>\\ a C a [ g a (s)ds\Y\ sup Y F a (d(x,z)). 

Jo zeY 



xeA n \A„ 



As m,n ^ oo, the above sum goes to zero. This proves that the sequence is Cauchy and hence 
convergent. The remaining claims follow as in Theorem 6.2.11 of pQ. □ 

3. Growth of Spatial Correlations 

The goal of this section is to prove Theorem l3.1l below which bounds the rate at which correlations 
can accumulate, under the influence of the dynamics, starting from a product state. 

3.1. The Main Result. Let be a normalized product state, i.e. 0, = &) xe j^£l x , where for 
each x, Q x is a state (not necessarily pure) for the systems at site x. We will denote by (•) the 
expectation with respect to f2, and prove 

Theorem 3.1. Let a > 0, G B a (A), and take 0, to be a normalized product state as described 
above. Given 1,7 C A with d(X,Y) > and local observables A G Ax and B G Ay, one has that 

(3.1) | (r t (AB)) - (r t (A)) (r t (B)) | < 4 \\B\\ \\F\\ ( \X\ + \Y\ ) G a (t) e ~ ad ^ Y \ 

Here 

(3-2) G a {t) = ° a + J FJ \Ma /" 9a(s)ds, 

JO 

and g a is the function which arises in the Lieb-Robinson estimate Theorem \2.1\ 

In the event that a = 0, the bound above does not decay. However, the estimate ()3.24j) below, 
which does decay, is valid. Moreover, a straight forward application of the techniques used below 
also provides estimates on the increase of correlations, due to the dynamics, for non-product states. 

We begin by writing the interaction $ as the sum of two terms, one of which decouples the 
interactions between observables supported near X and Y. 

3.1.1. Decoupling the Interaction: Consider two separated local observables, i.e., A G Ax and 
B eA Y with d(X, Y) > 0. Let 

(3.3) S AB := ft/eA: d(y,X) < d[XA ' ] 



2 

denote the ball centered at X with distance d(X,Y)/2 from Y. For any $ G B a (A), write 
(3.4) $ = $(1 - X a,b) + $XA,B =■ $i + $2, 

where for any Z C A 

1 if Z n S a ,b ^0andZn^ B ^0, 



(3-5) XA,b{ Z ) 
In this case, one has 



otherwise. 



Lemma 3.2. Let a > ; $ G B a {A), and consider any two separated local observables A G Ax and 
B G Ay with d(X, Y) > 0. Writing <3? = $i + <3? 2 , as in one may show that 

r\t\ r i 

(3-6) / [H 2 ,tW(0)\ ds < 2\\0\\G a (t) Yl E Fa(d(x,o)), 

oesupp(O) rtA: 

2d(x,o)>d(X,Y) 
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is valid for observables O G {A,B}. One may take 



(3.7) 



G a (t) 



C a + \\F a \ 



C a 



l$l 



9a(s) ds, 



where g a is the function from Theorem \2.1\ 
Proof. For O € {A B} and s > 0, 



(3.8) 



H 2 , r«(0) 



< 



E 



*(Z),r«(0) 



znS A , B ^9,znS c A s ^0 



as is clear from the definition of xa,b] see (|3-5(l . Applying Theorem 12.11 to each term above, we 
find that 

2g a {s)\\0\\mZ)\\ 



(3.9) 



*(Z),rW(0) 



< 



C„ 



£ E 

zgz oesupp(o) 



One may estimate the sums which appear above as follows: 



E E 

ZCA: ZgZ 



E 



\ 



E + E 

zeZ: zeZ: 



(3.10) 



Zr\S AtB jtQ,Zr\S% iB jtQ \z&Sa]b z€S% b J 

* E E E+ E E E. 

zGSa.b ^GS^ b Z^z,x zGS a s ^GSa.s Z3z,x 



and therefore, we have the bound 
(3.11) /"" [^ 2 ,r«(0) 

where 
(3.12) 



(is < 



gjjQj 



(Si + 5 2 ) / g a (s)ds 



5 i = E EE n*( z )ii E ^(^°)) 

2GSa,b x€S a b Z3z,x oGSUpp(O) 



and 
(3.13) 



5 * = E EE n*( z )ii E ^«(^°))- 

ze5 A,s xeS, A,s zbz,x oesupp(O) 
In the event that the observable O = A, then one may bound S\ by 

(3-14) S 1 < ||*|| a £ E Fa(d(z,x))J2 F a(d(z,y)) 

zeS A>B xes% B y ex 

< c a \m a E F «( d ( x 'y)) 

ieS^ B y£X 

and similarly, 

(3-15) S 2 < ||*|| a Y E Fa(d(z,x))Y F a(d(z,y)) 

z€S a,b x ^ s a.b y&x 

< WFaHHa E 52 F *Wz,y)) 

zdS A B yex 

An analogous bound holds in the case that O = B. We have proven (|3,6|) . □ 
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3.1.2. Proof of Theorem \S.l\ To prove Theorem 13.11 we will first provide an estimate which mea- 
sures the effect on the dynamics resulting from dropping certain interaction terms. 

Lemma 3.3. Let <3?o = 3>i + $2 be an interaction on A for which each of the dynamics {r/ }teR> 
for i £ {0, 1, 2}, exists as a strongly continuous group of *- automorphisms on A. Let {At}t£R be a 
differentiate family of quasi-local observables on A. The estimate 

,\t\ 



(3.16) \\Tp\At) -r^(A t )\\ < [ [H 2 ,rP(A s )] 

Jo 



+ 



f\d s A s )-r^(d s A s ) 



ds, 



holds for all ieR. Here, for each j 6 {0, 1, 2}, we denote by Hi the Hamiltonian corresponding to 



-4 by 



f{t) := rf\A t ) - ri l \A t ). 



Proof. Define the function / : 
(3.17) 

A simple calculation shows that / satisfies the following differential equation: 

(3.18) f{t) = i[H ,fit)] +l [H 2 ,r^iA t )] +r^id t A t )-ri 1 \d t A t ), 

subject to the boundary condition /(0) = 0. The first term appearing on the right hand side of 
(|3.18|) above is norm preserving, and therefore, Lemma |A. II implies that 

-1*1 



(3.19) 

as claimed. □ 



/(*) II < 



+ 



f\d s A s )-rPid s A s ) 



ds, 



We will now prove Theorem l3.ll Denote by Bt := B — (r^(i?)), and observe that proving ()3.1I) is 
equivalent to bounding |(Tt(AE^))|. Write $ = $1 + $2, as is done in ()3.4j) . One easily sees that 
$1 decouples A from B, i.e., 

(3.20) (r^iAB)) = (r t (1) (A) } (r t (1) (i?) }. 

Here, again, we have denoted by ri the time evolution corresponding to $1. It is clear that 



(3.21) 



\(n(AB t ))\ < (r^iABt)} + (r t (AB t ) - rr(ABt)) 



< 



\A\ 



r t iB)-T^\B 



+ 



Tt(ABt) — T^\AB- 



Moreover, the second term on the right hand side above can be further estimated by 



(3.22) 



T t iAB t )-T t (1) iAB 



< 2\\B\ 



T t iA)-T?\A) 



+ \\A\ 



rtiB t )-r^iB t ) 



Applying Lemma l3~3*l to the bounds we have found in (|3.21j) and (|3.22|) yields 



(3.23) 



{T t (AB t ))\ < 2\\A\ 



Hi,rPiB) 



ds + 21151 



1*1 



H 2 ,rPiA) 



ds. 



In fact, we are only using (|3.16|) in trivial situations where the second term, i.e., T s id s A s )—Ts 1 \d s A s 
is identically zero. Finally, using Lemma 13.21 we find an upper bound on \(rtiABt))\ of the form 



/ 



(3.24) 4 \\A\\ \\B\\ G a {t) 



\ 



E E Faidix,y)) + E E Faidix,y)) 



2d(x,y)>d(X,Y) 



y&Y ieA: 

2d(x,y)>d(X,Y) 



J 



Theorem 13. II readily follows from (|3,24|) above. 
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Appendix A. 

In this appendix, we recall a basic lemma about the growth of the solutions of first order, 
inhomogeneous differential equations. 

Let B be a Banach space. For each t £ M, let A(t) : B — > B be a linear operator, and denote by 
X(t) the solution of the differential equation 

(A.l) d t X{t) = A{t)X{t) 

with boundary condition X(0) = xq € B.We say that the family of operators A(t) is norm-preserving 
if for every xq € £>, the mapping j t : B — > B which associates xq — > X(t), i.e., Jt(xo) = X(t), satisfies 



(A.2) || 7t(x ) || = || x Q (I for all t G R. 

Some obvious examples are the case where B is a Hilbert space and A(t) is anti-hermitian for 
each t, or when B is an *-algebra of operators on a Hilbert space with a spectral norm and, for 
each t, A(t) is a derivation commuting with the *-operation. 

Lemma A.l. Let A(t), for t € R, be a family of norm preserving opeartors in some Banach space 
B. For any function B : R — > B, the solution of 

(A.3) 3 t Y{t) = A{t)Y{t) + B(t), 

with boundary condition Y(0) = yo, satisfies the bound 

(A.4) \\Y(t) - lt (y )\\ < I \\B{t')\\dt'. 

Jo 

Proof. For any t G R, let X(t) be the solution of 

(A.5) 8 t X(t) = A(t)X(t) 

with boundary condition X(0) = Xq, and let j t be the linear mapping which takes Xq to X(t). By 
variation of constants, the solution of the inhomogeneous equation (|A.3|) may be expressed as 

(A.6) Y(t) = 7 t(yo + ^(7s)- 1 (i?( S ))^). 

The estimate (|A.4|I follows from (|A.6|I as A(t) is norm preserving. □ 
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